Multi-layered genome defences in bacteria.

Bacteria have evolved a variety of defence mechanisms to protect against mobile genetic elements, including restriction-modification systems and CRISPR-Cas. In recent years, dozens of previously unknown defence systems (DSs) have been discovered. Notably, diverse DSs often coexist within the same genome, and some co-occur at frequencies significantly higher than would be expected by chance, implying potential synergistic interactions. Recent studies have provided evidence of defence mechanisms that enhance or complement one another. Here, we review the interactions between DSs at the mechanistic, regulatory, ecological and evolutionary levels.

Pseudomonas aeruginosa's adaptive trajectory: diverse origins, convergent paths.

Does genetic background contribute to populations following the same or divergent adaptive trajectories? A recent study by Filipow et al. evolved multiple genetically distinct Pseudomonas aeruginosa strains to an artificial cystic fibrosis lung sputum media. The strains adapted at different rates but converged on similar phenotypes despite their initial diversity.

Mutation bias and adaptation in bacteria.

Genetic mutation, which provides the raw material for evolutionary adaptation, is largely a stochastic force. However, there is ample evidence showing that mutations can also exhibit strong biases, with some mutation types and certain genomic positions mutating more often than others. It is becoming increasingly clear that mutational bias can play a role in determining adaptive outcomes in bacteria in both the laboratory and the clinic. As such, understanding the causes and consequences of mutation bias can help microbiologists to anticipate and predict adaptive outcomes. In this review, we provide an overview of the mechanisms and features of the bacterial genome that cause mutational biases to occur. We then describe the environmental triggers that drive these mechanisms to be more potent and outline the adaptive scenarios where mutation bias can synergize with natural selection to define evolutionary outcomes. We conclude by describing how understanding mutagenic genomic features can help microbiologists predict areas sensitive to mutational bias, and finish by outlining future work that will help us achieve more accurate evolutionary forecasts.

Mutational hotspots lead to robust but suboptimal adaptive outcomes in certain environments.

The observed mutational spectrum of adaptive outcomes can be constrained by many factors. For example, mutational biases can narrow the observed spectrum by increasing the rate of mutation at isolated sites in the genome. In contrast, complex environments can shift the observed spectrum by defining fitness consequences of mutational routes. We investigate the impact of different nutrient environments on the evolution of motility in Pseudomonas fluorescens Pf0-2x (an engineered non-motile derivative of Pf0-1) in the presence and absence of a strong mutational hotspot. Previous work has shown that this mutational hotspot can be built and broken via six silent mutations, which provide rapid access to a mutation that rescues swimming motility and confers the strongest swimming phenotype in specific environments. Here, we evolved a hotspot and non-hotspot variant strain of Pf0-2x for motility under nutrient-rich (LB) and nutrient-limiting (M9) environmental conditions. We observed the hotspot strain consistently evolved faster across all environmental conditions and its mutational spectrum was robust to environmental differences. However, the non-hotspot strain had a distinct mutational spectrum that changed depending on the nutrient environment. Interestingly, while alternative adaptive mutations in nutrient-rich environments were equal to, or less effective than, the hotspot mutation, the majority of these mutations in nutrient-limited conditions produced superior swimmers. Our competition experiments mirrored these findings, underscoring the role of environment in defining both the mutational spectrum and the associated phenotype strength. This indicates that while mutational hotspots working in concert with natural selection can speed up access to robust adaptive mutations (which can provide a competitive advantage in evolving populations), they can limit exploration of the mutational landscape, restricting access to potentially stronger phenotypes in specific environments.

Evolutionary innovation through transcription factor rewiring in microbes is shaped by levels of transcription factor activity, expression, and existing connectivity.

The survival of a population during environmental shifts depends on whether the rate of phenotypic adaptation keeps up with the rate of changing conditions. A common way to achieve this is via change to gene regulatory network (GRN) connections-known as rewiring-that facilitate novel interactions and innovation of transcription factors. To understand the success of rapidly adapting organisms, we therefore need to determine the rules that create and constrain opportunities for GRN rewiring. Here, using an experimental microbial model system with the soil bacterium Pseudomonas fluorescens, we reveal a hierarchy among transcription factors that are rewired to rescue lost function, with alternative rewiring pathways only unmasked after the preferred pathway is eliminated. We identify 3 key properties-high activation, high expression, and preexisting low-level affinity for novel target genes-that facilitate transcription factor innovation. Ease of acquiring these properties is constrained by preexisting GRN architecture, which was overcome in our experimental system by both targeted and global network alterations. This work reveals the key properties that determine transcription factor evolvability, and as such, the evolution of GRNs.

IκBζ is an essential mediator of immunity to oropharyngeal candidiasis.

Fungal infections are a global threat; yet, there are no licensed vaccines to any fungal pathogens. Th17 cells mediate immunity to Candida albicans, particularly oropharyngeal candidiasis (OPC), but essential downstream mechanisms remain unclear. In the murine model of OPC, IκBζ (Nfkbiz, a non-canonical NF-κB transcription factor) was upregulated in an interleukin (IL)-17-dependent manner and was essential to prevent candidiasis. Deletion of Nfkbiz rendered mice highly susceptible to OPC. IκBζ was dispensable in hematopoietic cells and acted partially in the suprabasal oral epithelium to control OPC. One prominent IκBζ-dependent gene target was ß-defensin 3 (BD3) (Defb3), an essential antimicrobial peptide. Human oral epithelial cells required IκBζ for IL-17-mediated induction of BD2 (DEFB4A, human ortholog of mouse Defb3) through binding to the DEFB4A promoter. Unexpectedly, IκBζ regulated the transcription factor Egr3, which was essential for C. albicans induction of BD2/DEFB4A. Accordingly, IκBζ and Egr3 comprise an antifungal signaling hub mediating mucosal defense against oral candidiasis.

SARS-CoV-2 ORF8 Mediates Signals in Macrophages and Monocytes through MyD88 Independently of the IL-17 Receptor.

SARS-CoV-2 has caused an estimated 7 million deaths worldwide to date. A secreted SARS-CoV-2 accessory protein, known as open reading frame 8 (ORF8), elicits inflammatory pulmonary cytokine responses and is associated with disease severity in COVID-19 patients. Recent reports proposed that ORF8 mediates downstream signals in macrophages and monocytes through the IL-17 receptor complex (IL-17RA, IL-17RC). However, generally IL-17 signals are found to be restricted to the nonhematopoietic compartment, thought to be due to rate-limiting expression of IL-17RC. Accordingly, we revisited the capacity of IL-17 and ORF8 to induce cytokine gene expression in mouse and human macrophages and monocytes. In SARS-CoV-2-infected human and mouse lungs, IL17RC mRNA was undetectable in monocyte/macrophage populations. In cultured mouse and human monocytes and macrophages, ORF8 but not IL-17 led to elevated expression of target cytokines. ORF8-induced signaling was fully preserved in the presence of anti-IL-17RA/RC neutralizing Abs and in Il17ra-/- cells. ORF8 signaling was also operative in Il1r1-/- bone marrow-derived macrophages. However, the TLR/IL-1R family adaptor MyD88, which is dispensable for IL-17R signaling, was required for ORF8 activity yet MyD88 is not required for IL-17 signaling. Thus, we conclude that ORF8 transduces inflammatory signaling in monocytes and macrophages via MyD88 independently of the IL-17R.

Transient mutation bias increases the predictability of evolution on an empirical genotype-phenotype landscape.

Predicting how a population will likely navigate a genotype-phenotype landscape requires consideration of selection in combination with mutation bias, which can skew the likelihood of following a particular trajectory. Strong and persistent directional selection can drive populations to ascend toward a peak. However, with a greater number of peaks and more routes to reach them, adaptation inevitably becomes less predictable. Transient mutation bias, which operates only on one mutational step, can influence landscape navigability by biasing the mutational trajectory early in the adaptive walk. This sets an evolving population upon a particular path, constraining the number of accessible routes and making certain peaks and routes more likely to be realized than others. In this work, we employ a model system to investigate whether such transient mutation bias can reliably and predictably place populations on a mutational trajectory to the strongest selective phenotype or usher populations to realize inferior phenotypic outcomes. For this we use motile mutants evolved from ancestrally non-motile variants of the microbe Pseudomonas fluorescens SBW25, of which one trajectory exhibits significant mutation bias. Using this system, we elucidate an empirical genotype-phenotype landscape, where the hill-climbing process represents increasing strength of the motility phenotype, to reveal that transient mutation bias can facilitate rapid and predictable ascension to the strongest observed phenotype in place of equivalent and inferior trajectories. This article is part of the theme issue 'Interdisciplinary approaches to predicting evolutionary biology'.

Causes of Death Among Individuals With Systemic Lupus Erythematosus by Race and Ethnicity: A Population-Based Study.

OBJECTIVE: Non-White populations are at higher risk of developing systemic lupus erythematosus (SLE) and have more severe outcomes, including mortality. The present study was undertaken to examine how specific causes of death vary by race and ethnicity, including Asian and Hispanic individuals. METHODS: The California Lupus Surveillance Project included SLE cases identified among residents of San Francisco County, CA during January 1, 2007 to December 31, 2009. Cases were matched to the National Death Index over a 10-year period. Logistic regression examined age-adjusted differences in causes of death by race, ethnicity, and sex. Age-standardized mortality ratios between individuals with SLE and the corresponding general population were calculated for the leading cause of death, and observed versus expected deaths were estimated. RESULTS: The study included 812 individuals of White (38%), Asian (36%), Black (20%), and mixed/other/unknown (5%) race; 15% identified as Hispanic. One hundred thirty-five deaths were recorded, with a mean ± SD age at death of 62.2 ± 15.6 years. Cardiovascular disease (CVD) was the leading cause of death overall (33%), and across all racial and ethnic groups, followed by rheumatic disease (18%) and hematologic/oncologic conditions (18%). CVD as the underlying cause of death was 3.63 times higher among SLE cases than in the general population. CVD deaths for those with SLE were nearly 4 and 6 times higher for Asian and Hispanic individuals with SLE, respectively, compared to the general population. CONCLUSION: Individuals with SLE experience a disproportionate burden of CVD mortality compared to the general population, which is magnified for Asian and Hispanic groups.

A gut-oral microbiome-driven axis controls oropharyngeal candidiasis through retinoic acid.

A side effect of antibiotics is outgrowth of the opportunistic fungus Candida albicans in the oropharynx (oropharyngeal candidiasis, OPC). IL-17 signaling is vital for immunity to OPC, but how the microbiome impacts antifungal immunity is not well understood. Mice in standard specific pathogen-free (SPF) conditions are resistant to OPC, whereas we show that germ-free (GF) or antibiotic-treated mice are susceptible. Oral type 17 cells and IL-17-dependent responses were impaired in antibiotic-treated and GF mice. Susceptibility could be rescued in GF mice by mono-colonization with segmented filamentous bacterium (SFB), an intestine-specific constituent of the microbiota. SFB protection was accompanied by restoration of oral IL-17+CD4+ T cells and gene signatures characteristic of IL-17 signaling. Additionally, RNA-Seq revealed induction of genes in the retinoic acid (RA) and RA receptor-α (RARα) pathway. Administration of RA rescued immunity to OPC in microbiome-depleted or GF mice, while RAR inhibition caused susceptibility in immunocompetent animals. Surprisingly, immunity to OPC was independent of serum amyloids. Moreover, RAR inhibition did not alter oral type 17 cytokine levels. Thus, mono-colonization with a component of the intestinal microflora confers protection against OPC by type 17 and RA/RARα, which act in parallel to promote antifungal immunity. In principle, manipulation of the microbiome could be harnessed to maintain antifungal immunity.

Environmental and societal factors associated with COVID-19-related death in people with rheumatic disease: an observational study.

Background: Differences in the distribution of individual-level clinical risk factors across regions do not fully explain the observed global disparities in COVID-19 outcomes. We aimed to investigate the associations between environmental and societal factors and country-level variations in mortality attributed to COVID-19 among people with rheumatic disease globally. Methods: In this observational study, we derived individual-level data on adults (aged 18-99 years) with rheumatic disease and a confirmed status of their highest COVID-19 severity level from the COVID-19 Global Rheumatology Alliance (GRA) registry, collected between March 12, 2020, and Aug 27, 2021. Environmental and societal factors were obtained from publicly available sources. The primary endpoint was mortality attributed to COVID-19. We used a multivariable logistic regression to evaluate independent associations between environmental and societal factors and death, after controlling for individual-level risk factors. We used a series of nested mixed-effects models to establish whether environmental and societal factors sufficiently explained country-level variations in death. Findings: 14 044 patients from 23 countries were included in the analyses. 10 178 (72·5%) individuals were female and 3866 (27·5%) were male, with a mean age of 54·4 years (SD 15·6). Air pollution (odds ratio 1·10 per 10 µg/m3 [95% CI 1·01-1·17]; p=0·0105), proportion of the population aged 65 years or older (1·19 per 1% increase [1·10-1·30]; p<0·0001), and population mobility (1·03 per 1% increase in number of visits to grocery and pharmacy stores [1·02-1·05]; p<0·0001 and 1·02 per 1% increase in number of visits to workplaces [1·00-1·03]; p=0·032) were independently associated with higher odds of mortality. Number of hospital beds (0·94 per 1-unit increase per 1000 people [0·88-1·00]; p=0·046), human development index (0·65 per 0·1-unit increase [0·44-0·96]; p=0·032), government response stringency (0·83 per 10-unit increase in containment index [0·74-0·93]; p=0·0018), as well as follow-up time (0·78 per month [0·69-0·88]; p<0·0001) were independently associated with lower odds of mortality. These factors sufficiently explained country-level variations in death attributable to COVID-19 (intraclass correlation coefficient 1·2% [0·1-9·5]; p=0·14). Interpretation: Our findings highlight the importance of environmental and societal factors as potential explanations of the observed regional disparities in COVID-19 outcomes among people with rheumatic disease and lay foundation for a new research agenda to address these disparities. Funding: American College of Rheumatology and European Alliance of Associations for Rheumatology.

Arid5a Mediates an IL-17-Dependent Pathway That Drives Autoimmunity but Not Antifungal Host Defense.

IL-17 contributes to the pathogenesis of certain autoimmune diseases, but conversely is essential for host defense against fungi. Ab-based biologic drugs that neutralize IL-17 are effective in autoimmunity but can be accompanied by adverse side effects. Candida albicans is a commensal fungus that is the primary causative agent of oropharyngeal and disseminated candidiasis. Defects in IL-17 signaling cause susceptibility to candidiasis in mice and humans. A key facet of IL-17 receptor signaling involves RNA-binding proteins, which orchestrate the fate of target mRNA transcripts. In tissue culture models we showed that the RNA-binding protein AT-rich interaction domain 5A (Arid5a) promotes the stability and/or translation of multiple IL-17-dependent mRNAs. Moreover, during oropharyngeal candidiasis, Arid5a is elevated within the oral mucosa in an IL-17-dependent manner. However, the contribution of Arid5a to IL-17-driven events in vivo is poorly defined. In this study, we used CRISPR-Cas9 to generate mice lacking Arid5a. Arid5a -/- mice were fully resistant to experimental autoimmune encephalomyelitis, an autoimmune setting in which IL-17 signaling drives pathology. Surprisingly, Arid5a -/- mice were resistant to oropharyngeal candidiasis and systemic candidiasis, similar to immunocompetent wild-type mice and contrasting with mice defective in IL-17 signaling. Therefore, Arid5a-dependent signals mediate pathology in autoimmunity and yet are not required for immunity to candidiasis, indicating that selective targeting of IL-17 signaling pathway components may be a viable strategy for development of therapeutics that spare IL-17-driven host defense.

Natural selection on crosstalk between gene regulatory networks facilitates bacterial adaptation to novel environments.

At the level of the gene, mutation is the raw material for natural selection. However, at the level of the gene regulatory network (GRN), variation is revealed to selection via promiscuous regulator activity ('crosstalk'), which creates opportunities for genetic innovation that can facilitate adaptation. Many genetic and environmental features can contribute to increasing potential for crosstalk by facilitating non-cognate interactions between regulatory elements. If a novel interaction provides a fitness benefit, rewired GRNs with strengthened affinity for newly forged connections can be selected. Here, we identify factors that facilitate opportunities for crosstalk and rewiring between GRNs, consider whether features of some GRNs make them more 'rewireable' than others and if these features might constrain evolution towards convergent outcomes. We explore patterns from laboratory and natural microbial populations that show changes within GRNs during adaptation. Finally, we discuss the prospects and open questions in the field.

Characteristics associated with poor COVID-19 outcomes in individuals with systemic lupus erythematosus: data from the COVID-19 Global Rheumatology Alliance.

AIM: To determine characteristics associated with more severe outcomes in a global registry of people with systemic lupus erythematosus (SLE) and COVID-19. METHODS: People with SLE and COVID-19 reported in the COVID-19 Global Rheumatology Alliance registry from March 2020 to June 2021 were included. The ordinal outcome was defined as: (1) not hospitalised, (2) hospitalised with no oxygenation, (3) hospitalised with any ventilation or oxygenation and (4) death. A multivariable ordinal logistic regression model was constructed to assess the relationship between COVID-19 severity and demographic characteristics, comorbidities, medications and disease activity. RESULTS: A total of 1606 people with SLE were included. In the multivariable model, older age (OR 1.03, 95% CI 1.02 to 1.04), male sex (1.50, 1.01 to 2.23), prednisone dose (1-5 mg/day 1.86, 1.20 to 2.66, 6-9 mg/day 2.47, 1.24 to 4.86 and ≥10 mg/day 1.95, 1.27 to 2.99), no current treatment (1.80, 1.17 to 2.75), comorbidities (eg, kidney disease 3.51, 2.42 to 5.09, cardiovascular disease/hypertension 1.69, 1.25 to 2.29) and moderate or high SLE disease activity (vs remission; 1.61, 1.02 to 2.54 and 3.94, 2.11 to 7.34, respectively) were associated with more severe outcomes. In age-adjusted and sex-adjusted models, mycophenolate, rituximab and cyclophosphamide were associated with worse outcomes compared with hydroxychloroquine; outcomes were more favourable with methotrexate and belimumab. CONCLUSIONS: More severe COVID-19 outcomes in individuals with SLE are largely driven by demographic factors, comorbidities and untreated or active SLE. Patients using glucocorticoids also experienced more severe outcomes.

A mutational hotspot that determines highly repeatable evolution can be built and broken by silent genetic changes.

Mutational hotspots can determine evolutionary outcomes and make evolution repeatable. Hotspots are products of multiple evolutionary forces including mutation rate heterogeneity, but this variable is often hard to identify. In this work, we reveal that a near-deterministic genetic hotspot can be built and broken by a handful of silent mutations. We observe this when studying homologous immotile variants of the bacteria Pseudomonas fluorescens, AR2 and Pf0-2x. AR2 resurrects motility through highly repeatable de novo mutation of the same nucleotide in >95% lines in minimal media (ntrB A289C). Pf0-2x, however, evolves via a number of mutations meaning the two strains diverge significantly during adaptation. We determine that this evolutionary disparity is owed to just 6 synonymous variations within the ntrB locus, which we demonstrate by swapping the sites and observing that we are able to both break (>95% to 0%) and build (0% to 80%) a deterministic mutational hotspot. Our work reveals a key role for silent genetic variation in determining adaptive outcomes.

The m6A reader IMP2 directs autoimmune inflammation through an IL-17- and TNFα-dependent C/EBP transcription factor axis.

Excessive cytokine activity underlies many autoimmune conditions, particularly through the interleukin-17 (IL-17) and tumor necrosis factor-α (TNFα) signaling axis. Both cytokines activate nuclear factor κB, but appropriate induction of downstream effector genes requires coordinated activation of other transcription factors, notably, CCAAT/enhancer binding proteins (C/EBPs). Here, we demonstrate the unexpected involvement of a posttranscriptional "epitranscriptomic" mRNA modification [N6-methyladenosine (m6A)] in regulating C/EBPß and C/EBPδ in response to IL-17A, as well as IL-17F and TNFα. Prompted by the observation that C/EBPß/δ-encoding transcripts contain m6A consensus sites, we show that Cebpd and Cebpb mRNAs are subject to m6A modification. Induction of C/EBPs is enhanced by an m6A methylase "writer" and suppressed by a demethylase "eraser." The only m6A "reader" found to be involved in this pathway was IGF2BP2 (IMP2), and IMP2 occupancy of Cebpd and Cebpb mRNA was enhanced by m6A modification. IMP2 facilitated IL-17-mediated Cebpd mRNA stabilization and promoted translation of C/EBPß/δ in response to IL-17A, IL-17F, and TNFα. RNA sequencing revealed transcriptome-wide IL-17-induced transcripts that are IMP2 influenced, and RNA immunoprecipitation sequencing identified the subset of mRNAs that are directly occupied by IMP2, which included Cebpb and Cebpd Lipocalin-2 (Lcn2), a hallmark of autoimmune kidney injury, was strongly dependent on IL-17, IMP2, and C/EBPß/δ. Imp2-/- mice were resistant to autoantibody-induced glomerulonephritis (AGN), showing impaired renal expression of C/EBPs and Lcn2 Moreover, IMP2 deletion initiated only after AGN onset ameliorated disease. Thus, posttranscriptional regulation of C/EBPs through m6A/IMP2 represents a previously unidentified paradigm of cytokine-driven autoimmune inflammation.

The racial/ethnic and sociocultural aspects of the pandemic in rheumatology.

The disproportionate impact of coronavirus-2019 (COVID-19) on communities of color is gaining global attention. Current research demonstrates that historically marginalized populations are experiencing disproportionate levels of SARS-Cov-2 infection and adverse clinical outcomes. However, research examining whether COVID-19 outcomes vary by race and ethnicity within the rheumatic disease population is limited. This paper will review data showing how SARS-CoV-2 infection has differentially affected racial and ethnic minorities in the general population and those with rheumatic disease. We will also highlight disparities in rheumatic disease risk and severity that existed prior to the pandemic, and discuss recent work examining severe outcomes of COVID-19 in rheumatic disease patients by race and ethnicity. Finally, we propose several actionable steps for the rheumatology community to address COVID-19 health disparities, which may have long-term effects on patients with rheumatic disease.

Association of Race and Ethnicity With COVID-19 Outcomes in Rheumatic Disease: Data From the COVID-19 Global Rheumatology Alliance Physician Registry.

OBJECTIVE: Racial/ethnic minorities experience more severe outcomes of coronavirus disease 2019 (COVID-19) in the general US population. This study was undertaken to examine the association between race/ethnicity and COVID-19 hospitalization, ventilation status, and mortality in people with rheumatic disease. METHODS: US patients with rheumatic disease and COVID-19 were entered into the COVID-19 Global Rheumatology Alliance physician registry between March 24, 2020 and August 26, 2020 were included. Race/ethnicity was defined as White, African American, Latinx, Asian, or other/mixed race. Outcome measures included hospitalization, requirement for ventilatory support, and death. Multivariable regression models were used to estimate odds ratios (ORs) and 95% confidence intervals (95% CIs) adjusted for age, sex, smoking status, rheumatic disease diagnosis, comorbidities, medication use prior to infection, and rheumatic disease activity. RESULTS: A total of 1,324 patients were included, of whom 36% were hospitalized and 6% died; 26% of hospitalized patients required mechanical ventilation. In multivariable models, African American patients (OR 2.74 [95% CI 1.90-3.95]), Latinx patients (OR 1.71 [95% CI 1.18-2.49]), and Asian patients (OR 2.69 [95% CI 1.16-6.24]) had higher odds of hospitalization compared to White patients. Latinx patients also had 3-fold increased odds of requiring ventilatory support (OR 3.25 [95% CI 1.75-6.05]). No differences in mortality based on race/ethnicity were found, though power to detect associations may have been limited. CONCLUSION: Similar to findings in the general US population, racial/ethnic minorities with rheumatic disease and COVID-19 had increased odds of hospitalization and ventilatory support. These results illustrate significant health disparities related to COVID-19 in people with rheumatic diseases. The rheumatology community should proactively address the needs of patients currently experiencing inequitable health outcomes during the pandemic.

Cancer cell lines show high heritability for motility but not generation time.

Tumour evolution depends on heritable differences between cells in traits affecting cell survival or replication. It is well established that cancer cells are genetically and phenotypically heterogeneous; however, the extent to which this phenotypic variation is heritable is far less well explored. Here, we estimate the broad-sense heritability (H 2) of two cell traits related to cancer hallmarks--cell motility and generation time--within populations of four cancer cell lines in vitro and find that motility is strongly heritable. This heritability is stable across multiple cell generations, with heritability values at the high end of those measured for a range of traits in natural populations of animals or plants. These findings confirm a central assumption of cancer evolution, provide a first quantification of the evolvability of key traits in cancer cells and indicate that there is ample raw material for experimental evolution in cancer cell lines. Generation time, a trait directly affecting cell fitness, shows substantially lower values of heritability than cell speed, consistent with its having been under directional selection removing heritable variation.